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Random bond disorder in antiferromagnetic spin- 1/2 Heisenberg ladders is investigated using the 
Stochastic Series Expansion Quantum Monte Carlo method. We find that the effects of individual 
bond impurities vary strongly, depending on their position on legs and rungs. We initially focus 
on how the distribution of local bond energies depends on the impurity concentration. Then we 
study how the phase diagram of even-leg ladders is affected by random bond doping. We observe 
Bose glass phases in two regimes {h' < h ^ hd and h" < h < hc2) and a Bose-Einstein Condensate 
in-between. Their presence are discussed in the relation to the local bond energies. 



I. Introduction: The effects of doping on spin-gap 
compounds has attracted considerable interest; follow- 
ing the discovery of new families of materials whose low- 
energy magnetic properties can be described by antifer- 
romagnetic (AF) Heisenberg ladder models'^'. It is known 
that AF spin-1/2 Heisenberg ladders with an even num- 
ber of legs are characterized by a singlet ground state. 
This symmetry of the ground state breaks down beyond 
a sufficiently high applied magnetic field or beyond a crit- 
ical concentration of random dopants. In the former case, 
Bose-Einstein condensation (BEC) of spin triplets occurs 
in a magnetic field window, hci < h < hc2, i-e. in an 
intermediate field regime which separates the spin-gap 
phase {h < hd) and a fully polarized phase (h > /ic2)!^ 
In the latter case, one distinguishes two types of random- 
ness: site and bond disorder. Site disorder induces local 
moments, which can give rise to long range order via an 
order- by-disorder mechanism.'^''^ In contrast, the pres- 
ence of bond disorder can destroy magnetic long range 
order through quantum localization, leading to Bose glass 
formation , which has recently been observed in spin gap 
ladder compounds^ ^. Many interesting phenomena are 
induced by the simult aneou s presence of both magnetic 
field and randomness .'^^^ For example, in doped spin 
gap ladders one observes a very rich phase diagram, with 
a sequence of magnetic field controlled phases, inc luding 
superfluidity, BEC, Bose glass and full polarization.'^HDiSl 
Specifically, recent experiments on the doped compound 
IPA-Cu(Clo.95Bro.o5)3 imply the existence of another 
Bose glass at low applied magnetic fields, before the ap- 
pearance of BEC. While some of these results have been 
modeled quantitatively using a bond-disordered Heisen- 
berg Hamiltonian^'^ it still remains desirable to better 
understand the effects of random doping on local observ- 
ables. 

In this paper, we study how bond disorder affects the 
quantum phase diagram of even-leg antiferromagnetic 
Heisenberg ladders in the presence of an applied mag- 
netic field. Bond disorder occurs when dopant ions re- 
place the ions which act as bridges between the mag- 
netic ions. For instance, bond disorder can be intro- 
duced in IFA-CuCla by a partial substitution of non- 
magnetic Br~ for the likewise non-magnetic Cl^, affect- 
ing the bond angles in the Cu-halogen-halogen-Cu su- 
perexchange pathwayd^ESUl, The enhanced strength of 



magnetic interactions on the affected bonds leads to a 
Bose glass when a sufficiently strong magnetic field is 
applied. Making use of Quantum Monte Carlo simula- 
tions we can address the local bond energies, uniform and 
staggered magnetization at ultralow temperatures. As 
discussed below, we find that with increasing magnetic 
field, BEC and Bose glass phases appear preceding the 
critical magnetic fields {hci and hc2) of the pure system. 
We study the influence of bond impurities on the bond 
energies Ei, of its neighbors. This allows us to charac- 
terize the effects of bond impurities on the condensation 
and localization of triplons, which occur in BEC and Bose 
glass phase respectively. In this context, a "superfluid" 
phase corresponds to the delocalization of bosons, and 
it is characterized by a non-zero magnetization per- 
pendicular to the magnetic field. The Bose glass phase 
occurs due to the localization of triplons. It is character- 
ized by a finite slope of the uniform magnetization rn„ 
and a vanishing order parameter to^. 

II. Model: We examine the bond-disordered anti- 
ferromagnetic Heisenberg model on 2- and 4-leg ladders 
described by the Hamiltonian 

i/ = ^J,,S,S,+/i^5f, (1) 

where Si is the spin-1/2 operator at site i, h \s the ap- 
plied magnetic field, and J^- denotes the nearest-neighbor 
coupling between the spins on sites i and j, taking 
the values J = 1 with probability 1 — p and J' with 
probability p. Periodic boundary conditions are used 
along the leg direction. The simulations are performed 
using the Stochastic Series Expansion (SSE) Quantum 
Monte Carlo (QMC) method based on the directed loop 
algorithnJ^. In order to access the regime of very low 
temperatures, we use a /3— doubling schem^i®. This way, 
maximum inverse temperatures up to /? = 2048 are ob- 
tained. 

The staggered magnetization is calculated from 
the staggered structure factor, 

Si^]^Y.^-^y^'{SfS^), (2) 

using TO^ = ^/s]-Jn. And the local bond energies are 
measured by Eb = Jij{SiSj). The staggered and 



2 



B 15 
(X 

10 



20 
15 

of 10 





1 ' 1 ' 1 ' 






Q- o all bonds 








- — ■ leg bonds 








• — • rung bonds 










Q--0 all bonds 
- — ■ outer-leg bonds 
• — • inner-leg bonds 
* — * outer-rung bonds 
T — I center-rung bonds 



(a) 

3 -0.234 -0.390 
-0.507 

-0.418 -6.0328 -0.414 

(b) 

. -0.416 _-0.414. 



-0.397; :-0.059 -C 



-0.387 -0.355 -0.372 



(c) 



-0.389 
-0.322 



-0,267, 
■0,424 



■0.395 
-0.329 



■0.340 
■0,424 



-0.438 -0.038 



-0.421 -0.414 



.-0.371 ,-0.372 



(d) 



[ -0.384 1-0.377 j-0. 
' -0.320 y-0.319y 



-0.357* -0,355 A 
11-0.348; i-0. 048 J-0. 
-0.421 -0.419 



-0.385 _-0.361_ -0.374 




■0.390 -0.310 -0.3 



(f) 



-0.359 -0.360 



398 
.0.342 



0.435 [ 
.0.337 



-0.34in -0.340 
0.395 -0.432 



-0.360 -0.358 



FIG. 1: (color online) (a) Distribution of local bond ener- 
gies in pure (a) 2-leg and (b) 4- leg AF spin-1/2 Heisenberg 
ladders. For each case, there are different contributions due 
to rung bonds and leg bonds. Simulations are performed at 
temperature T = J/2048, and lattice length = 256. 



FIG. 2: (color online) Effects of local bond impurities, de- 
noted by dashed lines, on the neighboring bond energies in 2- 
and 4-leg ladders. Simulations are performed at temperature 
T = J/2048, lattice length Lx ~ 256, with the coupling of the 
bond impurity J' = J/ 5. 



uniform m„ magnetizations are averaged over at least 
700 impurity realizations. An ultra-low temperature T — 
1/2048 was chosen in all simulations so that the relevant 
thermodynamic observables reflect true zero-temperature 
behavior. The lengths of the 2-leg ladders used are — 
104,128,160,200,256,320, and 400. And the lengths of 
the 4-leg ladders used are = 96, 128, 128, 160, 256, 320, 
and 400. Their thermodynamic limit, i.e. infinite length, 
is extrapolated via finite-size scaling. 

III. Bond-doped Heisenberg ladders: Pure anti- 
ferromagnetic even-leg ladders display a resonant valence 
bond ground stat^ with a finite gap A to triplet excita- 
tions, where the n spins on the same rung preferentially 
form local singlets. The value of the spin gap decreases 
as the number of legs increases.'^ An interesting question 
that arises concerns how the magnitudes of the bond en- 
ergies relate to the observed spin gap. 

Fig. [T] shows the local bond energy distribution in pure 
2- and 4-leg ladders with Jij = J. The simulations 
are performed at a very low-temperature T = J/2048 
on ladders with linear size = 256, measuring the 
local bond energies. In pure 2-leg ladders, the distri- 
bution function of these local bond energies shows two 
peaks^'^, corresponding to two the types of bonds. The 
rung-bond energies give rise to the lower-energy peak at 
Et « — 0.454J, whereas the leg bonds cause the higher- 
energy peak at Eb « —0.350 J. This observation is consis- 
tent with the fact that pure 2-leg ladders display a rung- 
singlet-dominated ground state. The local bond energies 
of pure 4-leg ladders fall into four categories, correspond- 
ing to four types of bonds, located on outer rungs, outer 
legs, inner legs, and center rungs respectively. They are 
arranged in ascending order of local bond energies, given 
by Eb « -0.388 J, -0.370J, -0.326 J, and -0.309 J. It is 
noted that the center-rung bonds have the highest ener- 
gies. In other words, the four rung spins preferentially 
form effective pairs of weakly coupled singlets on each 



rung. 

What happens when the ladders are bond-doped? For 
sufficiently low bond dilution concentrations, the average 
distance between the individual bond impurities is large. 
To understand this regime, it is useful to study the local 
effects of isolated impurity bonds. Fig. [2] shows the effect 
of an isolated impurity on its neighbor bonds. Here, one 
bond is replaced by a weaker coupling J' = J/5, and its 
position is varied. In the case of 2-leg ladders, the im- 
purity bond can be located on either a rung or on a leg. 
In Fig. [2]^ a) , it is observed that a leg-bond impurity sig- 
nificantly reduces the magnitude of the bond energy on 
the opposite leg bond. Additionally, it strongly enhances 
the bond energies on its two neighbor rungs, indicating 
that this type of bond impurity strengthens the forma- 
tion of rung singlets in its immediate neighborhood. In 
contrast, when the bond impurity is placed on a rung, it 
reduces the magnitude of the bond energies of the two 
closest rungs, see Figj2|^b), but increases the magnitude 
of the bond energies on the neighboring legs. 

For the case of 4-leg ladders, the preferential formation 
of pairs of effectively weakly coupled singlets along the 
4-site rungs has significant consequences when bond dop- 
ing is introduced. As discussed above, an approximate 
effective description of the ground state in a pure 4-leg 
ladder would be a product state of two weakly coupled 2- 
leg ladders. When a bond impurity is introduced on a leg, 
it thus causes an effect similar to the 2-leg ladder case, 
see Fig. [2jc) and (e), i.e. the opposite leg bond energy 
on the 2-leg subsystem on which the impurity is placed 
is reduced while bond energies on neighboring rungs are 
enhanced. The other 2-leg subsystem is barely affected 
by the leg impurity. In contrast, if the impurity is intro- 
duced on an outer rung, it causes a reduction of the rung 
neighbor bond energies on the same 2-leg ladder subsys- 
tem, see Fig. [2]jd). It also enhances the bond energies on 
the neighboring leg bonds, which again is similar to the 
effect in Figj2Fb). Last but not least, when the impurity 
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FIG. 3: (color online) (a) Distribution of local bond energies 
in bond-disordered (a) 2-leg and (b) 4-leg ladders. Simula- 
tions are performed at temperature T = J/ 2048, and lattice 
length = 256. 



FIG. 4: (color online) Staggered and uniform magnetization 
of pure and bond-disordered even-leg ladders. Bond couplings 
are randomly replaced by couplings J' = J/5. The values of 
the staggered magnetization are averaged over at least 700 
realizations. 



is placed on a center rung, it enhances the bond energies 
on the two neighboring rungs. However, its influence 
compared with the previous case is much less drastic, see 
Fig. if). 

In randomly bond doped ladders, the distribution of lo- 
cal bond energies displays characteristics stemming from 
the effects discussed above for isolated impurity bonds. 
Fig. |3] shows the distribution of local bond energies of 
doped (a) 2-leg and (b) 4-leg ladders. Here, we show 
three typical examples with {z = 4.2%, J' = J/5), 
{z = 8.3%, J' = J/5), and (z = 8.3%, J' = 0). 

For the 2-leg ladder, it is observed that, besides the 
main peaks of the pure case, there are three additional 
features. For instance, let us consider the case (z = 
8.3%, J' = J/^)- Here, the additional peaks appear at 
energies « —0.51 J, — 0.22J, and —0.05 J. These features 
originate from affected bonds in the vicinity of the bond 
impurities. The lowest-energy peak is the contribution of 
singlets next to leg-bond impurities, the second lowest- 
energy peak comes from leg bonds opposite to the impu- 
rity, and the third peak stems from the bond-impurities 
themselves on the rung and leg bonds. There is actually 
another feature in-between the two main peaks, which 
stems from rung singlets next to the rung-bond impuri- 
ties. 

In 4-leg ladders with {z = 8.9%, J' = J/5), the feature 
at —0.42 J stems from the rung bonds shown in Figs.[2|^c), 
(e), and (f) and outer- leg bonds in Fig. [2jd). The peak 
at Eb ss — 0.27J is corresponds to center leg couplings 
in Fig. [2jc). The slightly lower energy peak at £'(, w 
— 0.28J is a contribution of center rung bonds next to 
rung-bond impurities, see Fig.[2][f). The peak at —0.03 J 
is the contribution of bond-impurities. Clearly, the local 
bond energy distribution is much richer compared with 
the 2-leg ladder case, but the individual features can still 
be explained using a local description. 

The strength of local bond energies has direct conse- 
quences for the phase diagram when the randomly bond 



doped ladders are exposed to an applied magnetic field. 
Fig. |4] shows the quantum (low-temperature) phase dia- 
gram of pure and doped 2 and 4-leg ladders in the pres- 
ence of a uniform magnetic field h. In the pure case, for 
sufficiently small fields, both the uniform m„ and stag- 
gered magnetization remain. They become finite be- 
yond a lower critical field, hd = 0.5 for the 2-leg ladder 
and hci — 0.16 for the 4-leg ladder. Spins in pure 2-leg 
ladders are fully polarized at an upper critical magnetic 
field hc2 = 3, and simil arly , the upper critical field for 4- 
leg ladders is hc2 = 3.8.123 Fig. [5]ja) provides a schematic 
picture to explain the observation of the phase diagram 
found from the QMC simulations. At zero field, quan- 
tum spin fluctuations in the pure 2-leg Heisenberg ladder 
destroy conventional magnetic order. The RVB ground 
state can be approximated by a product state of rung 
spins forming singlets, separated from the lowest triplet 
excitation by a minimum excitation energy A. This is a 
non-magnetic spin liquid state. The spin gap is reduced 
as the magnetic field increases. It closes when the mag- 
netic field reaches the lower critical value hd, at which 
the system experiences a quantum phase transition from 
the spin liquid phase to a BEC of magnons.l^This phase 
is characterized by non-zero uniform and staggered mag- 
netizations. Ultimately, at large applied magnetic fields 
all the spins are polarized along the z-direction, i.e. for 
h > hc2- At this upper critical magnetic field, the sys- 
tem transitions from a BEC to a polarized paramagnet, 
and there is no longer coherence between spins. The 
staggered magnetization vanishes and the uniform mag- 
netization saturates at m„ = 1/2. 

In the doped case, these phase diagrams are modified 
by the presence of bond disorder. The staggered magneti- 
zation of doped 2 and 4-leg ladders are shown for respec- 
tive doping concentrations z = 4.2% and 4.5%. The bond 
couplings J are randomly replaced by J' = J/5. This 
particular choice is motivated by recent experiments,!^ 
but the results shown here are generic. Interestingly, 
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FIG. 5: (color online) Schematic illustration of the response 
of the uniform and staggered magnetizations to an applied 
magnetic field. In the pure case (a), the system remains in 
an RVB spin fluid state up to a lower critical field . Then 
it undergoes a BEG transition into the regime < h < hc2 ■ 
All spins are fully polarized beyond the saturation field hc2- In 
the doped case (b), a field scan reveals the following sequence 
of phases: I: RVB spin liquid; II: Bose glass phase; III: BEG; 
IV: Bose glass phase; V: fully polarized phase. 



close to both critical fields h' < h ^ hd and h" < h < hc2 
(Fig.|4jb) and (c)) we observe a finite slope of the uniform 
magnetization but a vanishing staggered magnetization. 
These are the characteristics of Bose glass phases. We ob- 
serve the phase transition of these Bose glass phases to 
BEC at a magnetic field smaller than hd for 2-leg ladders 
and greater than hci for 4-leg ladders. Observations of 
disorder induced phases near the originally critical fields, 
recently reported in three-dimensional dimer systems'^ 
suggest that BEC and Bose glass phases appear near the 
higher critical fields. Our finding differs in that for the 
bond-doped 2-leg ladder the Bose glass phases occur at 
both critical fields, depending on the magnitude of bond 
disorder. However the existence of Bose glass phases at 
both critical magnetic fields has recently been observed 
in the compound IPA-Cu(Clo.95Bro.o5)3.'^Below, we dis- 
cuss mechanisms derived from our model study which 
offer explanations for the experimental observations. 

In Fig. Sjjb), we show the five phases which appear 
in the QMC data of doped 2-leg ladders. At sufficiently 
small fields the system stays in the spin liquid regime 
I. A finite magnetic field strength is required to over- 
come the lowest singlet-triplet gap of bond-impurity cou- 
plings. As discussed above, the impurity rung-bonds are 
weaker, and therefore singlets on these bonds break first, 
when an external magnetic field is applied to the sys- 
tem. Beyond this first critical field, the uniform mag- 
netization becomes finite with increasing magnetic field 
while the staggered magnetization stays zero. This mag- 
netic field is sufficiently strong to gradually polarize spins 
on the impurity rung-bonds, but not strong enough to 



polarize the spins in the bulk. During this phase, the 
field-induced triplons stay localized on the impurity-rung 
bonds. There is no coherence among them. The system 
therefore forms a Bose glass phase which is manifested 
in region II. At magnetic field, h « hd, singlets on the 
bulk rungs break, triplons form on these rung-bonds and 
they interact with each other. We note that for any small 
amount of dilution impurities in a 2-leg ladder separate 
rung-bonds into clusters. Rung couplings on these clus- 
ters are generally greater than that in the pure case. 
Hence the result it is required a magnetic filed higher 
than hd to polarize these rung-bonds, FigQc). However, 
the above argument can not be applied to 4-leg ladders 
because at small dilution rung-bonds are not separated 
into clusters. Therefore a magnetic field smaller than hd 
is sufficient to polarize these rung -bonds, Figgb). The 
system undergoes a BEC phase transition, regime III. Re- 
call that there are still basically un-renormalized singlets 
next to the leg-bond impurities untouched. Their bond 
energies are relatively large compared with the original 
bond energies, see Fig.[2][a). Although the magnetic field 
has fully polarized the original singlets in this regime, it 
is not sufficiently strong to break these impurity induced 
higher-energy singlets. The system therefore undergoes 
another Bose glass phase transition into regime IV. Fi- 
nally, above hc2, all spins finally polarized, regime V. 
This sequence is observed in the QMC data shown in 
Fig.i 

IV. Conclusions: To summarize, we have studied 
the effects of bond impurities in even-leg antiferromag- 
netic spin-1/2 Heisenberg ladders. We find that depend- 
ing on the location of the bond impurities, bond energies 
on neighboring bonds are either enhanced or reduced. 
In general, bond impurities enhance bond energies con- 
nected with them by the same spins, and reduce the op- 
posite couplings. This effect is related to the so-called 
edge disorder effect .fSS We identify various types of impu- 
rity induced bond energy shifts in 2- and 4-leg ladders. 
In light of these results, we demonstrate the emergence of 
BEC and Bose glass phase close to both critical fields of 
the pure system. These results can be used to explain re- 
cent observations of a disorder induced Bose glass phase 
in IPA-Cu(Clo.95Bro.o5)3.'^ 
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